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Abstract 

We present lattice results for the spectrum of mesons containing one heavy quark and 
of baryons containing one or two heavy quarks. The calculation is done in the quenched 
approximation using the NRQCD formalism for the heavy quark. We analyze the de- 
pendence of the mass splittings on both the heavy and the light quark masses. Me- 
son P-state fine structure and baryon hyperfine splittings are resolved for the first 
time. We fix the b quark mass using both Mb and M\ b , and our best estimate is 
mff s \m^ IS ) = 4.35(10)(+|)(20) GeV. The spectrum, obtained by interpolation to nib, 
is compared with the experimental data. 



* Present address: Spatial Technologies, Boulder, CO, USA. 
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1 Introduction 



The spectrum and decays of hadrons containing b quarks will be measured in precision ex- 
periments at the B factories. It is therefore important to calculate the spectrum expected 
from QCD, both as a test of the theory and to predict the masses of states not yet observed. 
This paper reports on results of a lattice calculation of the heavy-light spectrum using the 
non-relativistic formulation of QCD (NRQCD) for heavy quarks and the tadpole-improved 
clover action for light quarks. This approach allows us to have better control over discretization 
errors in both the heavy and the light quark sectors. 

Lattice QCD allows us to investigate the dependence of the meson and baryon mass split- 
tings on the heavy and light quark masses. For this purpose we simulate three values of light 
quark masses in the range 0.8m s — 1.3m s , and six values of heavy quark masses in the range 
3 — 20 GeV. The NRQCD formalism is ideally suited to study such a wide range of heavy 
quark masses at 1/a = 1.92 GeV, the lattice spacing we use. For the light quarks we use 
the tadpole- improved clover action which has discretization errors of 0(a s a) and these are 
expected to be small at this lattice spacing. These improvements make it possible to perform 
reliable comparisons with both the experimental b spectrum and expectations based on Heavy 
Quark Symmetry. 

The phenomenological interest in the decay rates of hadrons containing b quarks stems 
from the important role they play in the determination of the Cabibbo-Kobayashi-Maskawa 
matrix elements. Two quantities that are used as input in the analyses of experimental data 
are mff 6 ! (m^ IS ) and the decay constants js and f Bs . Here we shall present results for the b 
quark mass, while the calculation of decay constants has already been reported in a companion 
paper @. 

This paper is organized as follows. In Sec. ^| we briefly review the experimental situation 
and provide a justification for the NRQCD approach to heavy quarks. The parameters used 
in the simulations are given in Sec. |3|. Sec. § describes the determination of the b quark mass. 
Our results on the heavy-light meson spectrum are presented in Sec. | along with a discussion 
of the spin-independent and spin-dependent mass splittings. Baryons containing one heavy 
and two light quarks are discussed in Sec. ||. In Sec. |7| we give a brief description of our results 
on baryons containing two (degenerate) heavy quarks and one light quark. This is followed by 
a determination of HQET parameters in Sec. []. Finally, we summarize our main conclusions 
in Sec. 

2 Phenomenological Background 

The NRQCD approach for simulating b quarks is justified because the typical velocity of 
the heavy quark is small, v/c ~ 0(Aq C d/M) ~ 0.05 — 0.1. This is corroborated by the 
experimental observation that all splittings are much smaller than the masses, and the hadron 
masses are dominated by the heavy quark mass. Thus a very natural picture of the heavy- 
light system is a "hydrogen atom" composed of the light degrees of freedom bound in the 
background of an almost static color source. Within this model one can distinguish between 
spin- independent splittings in the spectrum dominated by radial and/or orbital excitations 
of the light quark, and spin-dependent ones dominated by the spin-flip energy of the heavy 
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Figure 1: Experimental spin- independent mass splittings for hadrons with one heavy quark 
(h = b, c, or s) as a function of the spin-averaged meson mass M = (M# + 3M#*)/4 where 
H denotes a generic heavy meson. Squares denote the B s — Bd and the D s — Da splitting. 
Pluses stand for the spin-averaged X — A splitting (we have used the DELPHI measurement 
of Eft ||). The splitting between the A and the spin-averaged S state meson is denoted by 
crosses. Bursts denote the spin-averaged P — S splitting. 
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Figure 2: Experimental spin-dependent mass splittings for hadrons with one heavy quark. M 
is defined in Fig. [I]. Squares denote the S state hyperfine splitting for B, B s , D, D s , and K 
mesons. Diamonds denote the splitting between P states with ji = 3/2. These are known 
only for D, D s , and K mesons. The £* — S splitting for baryons is denoted by bursts. (For 
a discussion of the possibility that some of the S c and baryons have been misidentified 
experimentally see Ref. M). 
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Heavy quark parameters 


aM° 
n 


1.6 

2 


2.0 
2 


2.7 

2 


4.0 
1 


7.0 
1 


10.0 
1 


Light quark parameters 


K 


0.1369 
0.423(7) 


0.1375 
0.362(11) 


0.13808 
0.298(4) 



Table 1: A summary of the heavy and light quark mass parameters, aM° and k, used in the 
simulation and the resulting mass of pions composed of degenerate light quarks. We also list 
the values of the stability parameter n used in the heavy quark evolution 0. 



quark. These two types of splittings have distinct behavior as a function of the heavy quark 
mass. Spin-independent splittings survive the infinite heavy quark mass limit whereas the 
spin-dependent ones do not. 

The experimental data plotted in Fig. |T] show that the spin-independent splittings are 
often insensitive to the mass of the heavy quark. In fact one finds in many cases that the 
insensitivity persists down to the strange quark mass. Spin-dependent splittings, on the other 
hand, are found to increase with the inverse heavy quark mass as shown in Fig. |2]. An analysis 
with a phenomenologically determined potential is in agreement with these results, however 
there is considerable uncertainty in how to model the light degrees of freedom (see H and 
references therein). Simulations of lattice QCD using a non-relativistic formulation for heavy 
quarks provide estimates without resort to modeling. 

The NRQCD formulation has been discussed in |l|, §]. It has been very successful in the 
study of heavy quarkonia [0], and we apply it to predict the heavy-light spectrum here. Results 
using alternate formulations, static heavy quarks or standard (Wilson or clover) discretization 
of the Dirac operator mostly extrapolated from the charm region, can be found in J7|, ^, £| TU 



11| and we shall compare against them at appropriate places. 



3 Simulation Parameters 

The statistical sample consists of the same 102 quenched configurations, at /3 — 6.0 with 
lattice size 16 3 x 48, as used in our study of decay constants ||. The NRQCD action, the 
evolution equation for calculating the heavy quark propagator, the method used for setting 
the lattice scale, and the fixing of light and strange quark masses are also the same. The list 
of quark masses used in our simulation are reproduced in Table [I], and the operators used to 
study the various states, are given in Table |[ 

We estimate that the significant sources of systematic errors in this calculation are finite 
volume, finite lattice spacing, quenching, uncertainties in determining a, fixing the strange 
quark mass and perturbative corrections. For a lattice size of ~ 1.6 fm, finite volume effects 
are not expected to be significant for the lower lying S state mesons. However, there are 
indications that the wave functions for P states and the baryons are more extended M and 
finite size effects in these states should therefore be larger. We cannot comment on this as 



4 



state 


operator 


1 C 
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u <j h 


1 r> 
Pi 






— A 7 


3 Po 






ua-Ah 


3 Pl 






u a x A h 




(T) 




— /A A \ 7 

W (CjAj + CTj/Xi) h 


3 P 2 


(E) 




u (crjAj — <JjAj) h 


A 


(s* = 


+ 1/2) 


u c d 


S 


(s z = 


+1/2) 


(u c a z d h-\ — \/2u c a + d hi)/y3 


S* 


(s z = 


+3/2) 


u c a + d h-\ 


s* 


(S z = 


+1/2) 


(\/2u c a z d hi + u c a + d /ij,)/\/3 



Table 2: The operators used to study the various states, h stands for the two-component heavy 
quark spinor, u and d for the upper two components of two flavors of light quark spinors. A is 
the ordinary derivative in the Coulomb gauge. The symbols and h± stand for the s z = +| 
and — | components of h respectively. The baryon operators for s z < are constructed from 
the corresponding s z > operators by interchanging a + «-> cr_ and The H baryons are 

obtained by replacing one of the light flavors in S by an s, and the f2 by replacing both light 
quarks by ss. For the heavy-heavy-light baryons, the operators are identical except u and d 
are to be interpreted as two flavors of heavy quarks and h as the light or s quark. The 3 P2 
states decompose, under the cubic group, into two representations labeled T and E. Our j — 2 
P states are spin-averaged over both lattice representations: 3 P 2 = [3 3 P 2 (T) + 2 3 P 2 (E)]/5. 



we have results on only one lattice volume. The 0(a s a) error associated with the tadpole 
improved clover light fermions is expected to be a few percent at this (3 |T2]|. A detailed study 



of the scaling behavior of the heavy-light spectrum is discussed in Ref. ||13|| . Quenching errors 
remain unknown. However, since the B spectrum is dominated by the light quark degrees of 
freedom, we expect that using light spectroscopic quantities to fix a compensates for part of 
this uncertainty. 

The central value of lattice scale we use is 1/a = 1.92(7) GeV as obtained from M p . To 
estimate the systematic error in this we repeat our bootstrap analyses with 1/a = 1.8 and 2 
GeV as discussed in [0. We obtain Ki = 0.13917(9), corresponding to the light quark mass 
mi = 1/2 (m u + md), by linearly extrapolating M^/M p to 137 2 / 770 2 . We cannot resolve a 
curvature in the light quark mass dependence, and do not assign a systematic error in k\. To 
determine the strange quark mass, we use three different methods. By fixing the ratio M|-/M 2 
to its physical value, we obtain k s = 0.13755(13). Using the ratios Mk*/M p and M^/Mp, gives 
k s = 0.13719(25) and 0.13717(25) respectively. Since the latter two agree within errors, we 
only give the results using Mk and Mk*- For our final results, we use n s from Mk, and 
determine the systematic error using k s from K*. 

In our final results, the first error we quote comes from a bootstrap analysis using a -1 = 
1.92(7) GeV, the second from the scale uncertainty, and where applicable, the third from the 
uncertainty in the strange quark mass. We comment on the uncertainty due to using 1-loop 
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perturbative expressions and in fixing the b quark mass below. 

A summary of some of the important features of the raw lattice data are as follows, (i) The 
data at aM° = 7.0 and 10.0 are not as reliable as that for aM° < 4.0 (there are no clear 
plateaux in the effective mass plots). They are, therefore, used only in the estimation of 
HQET parameters, where we have chosen states and operators with the best signal, (ii) The 
calculation of P state correlation functions has been done for only aM° = 1.6, 2.0, 2.7. 

Lastly, we fix the bare b quark mass aM" as follows. In NRQCD and the static theory, 
E s i m , the rate of exponential fall-off of the heavy-light meson correlators, is not the meson 
mass, but is related to it by the shift, 

A = Mhadron _ -^sim = Z m M° — Eq . (1) 

Here Z m is the renormalization constant connecting the bare quark mass to the pole mass, 
and Eq is the shift in the energy of the quark. As discussed in more detail in Ref. f2|, we 
employ three different methods to calculate the meson mass: (i) M^ n extracted directly from 
the dispersion relation of the heavy-light meson; (ii) M pert obtained by evaluating the mass 
shift A perturbatively; and (iii) M' using the A obtained from the dispersion relation of the 
heavy-heavy meson at the same aM° value. The perturbative results for Z m , Eq and A are 
given in Table |||.[] 



aM° 


n 


7 


aE 


aA 


1.6 


2 


1.18 


0.23 


1.64 


2.0 


2 


1.14 


0.28 


2.02 


2.7 


2 


1.09 


0.27 


2.68 


4.0 


1 


1.05 


0.27 


3.90 


7.0 


1 


1.00 


0.28 


6.74 


10.0 


1 


0.98 


0.28 


9.54 



Table 3: The stability parameter n and the 1-loop perturbative estimates of the mass renor- 
malization constant Z m , the zero point shift of the heavy quark energy Eq, and the mass shift 
A = Z m M° — E Q using the q* calculated with the Hornbostel-Lepage procedure fl4j]. Errors 



associated with numerical integration of the 1-loop expressions are insignificant compared to 
other systematic errors. 

In the perturbative analyses, we use a s = ap defined in Ref. ||15|| . The relevant scale 
q* at which to evaluate the running coupling ap is chosen separately for each process using 
an extension [14] of the Lepage-Mackenzie scale-setting prescription fl6|]. The choice of scale 



advocated in the original Lepage-Mackenzie scheme eliminates the 0(a 2 s ) correction in the 
bubble summation approximation. This procedure can fail, however, when the one-loop con- 
tribution becomes small. Hornbostel and Lepage [T4] have recently extended the method to 



overcome this difficulty by taking into account higher-order terms in the bubble summation 



lr The perturbative calculations have been done for a slightly different discretization of F^, i.e. a four 
leaf clover rather than the two leaf version used in the evolution equation. We expect the difference to be 
insignificant. 
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K 


s 




aM® 


a E ■ 


aM kin 


aM' 


aMp Crt 


n E ■ 

Uj - LJ sim 


aM kin 


aM' 


aMpcrt 


1.6 


0.427(7) 


2.16(13) 


2.08(3) 


2.07(4) 


0.474(4) 


2.21(9) 


2.13(3) 


2.11(2) 


2.0 


0.443(8) 


2.57(18) 


2.46(4) 


2.46(2) 


0.490(4) 


2.63(13) 


2.50(4) 


2.51(1) 


2.7 


0.459(7) 


3.30(30) 


3.15(7) 


3.14(2) 


0.504(4) 


3.35(21) 


3.20(7) 


3.18(1) 


4.0 


0.468(8) 


4.76(65) 


4.46(11) 


4.37(7) 


0.513(5) 


4.74(43) 


4.50(11) 


4.41(7) 


7.0 


0.469(9) 


8.9(24) 




7.21(21) 


0.516(5) 


8.5(15) 




7.26(21) 


10. 


0.471(9) 


15(7) 




10.01(35) 


0.515(5) 


13(4) 




10.1(3) 



Table 4: E S [ m and pseudoscalar meson masses in lattice units extrapolated/interpolated to K\ 
and k s . Meson masses have been calculated from the heavy-light dispersion relation (M^ Q ), 
using A from heavy-heavy spectroscopy (M ; ), and from perturbation theory (M pert ). 



approximation. Their extension reduces to the original Lepage-Mackenzie prescription when 
the one-loop term is not small due to large cancelling contributions. 

The perturbative series for Z m has an infra-red renormalon ambiguity |l7j, which is typi- 
cally characterized by an uncertainty of 0(Aqcd/M). Since this is comparable to the entire 
0(a s ) correction, we shall use the latter as the estimate of the perturbative error in the 
determination of M pole . 

All three methods for estimating the B meson mass give compatible results for aM® < 4 as 
shown in Table |j. These estimates differ slightly from those in Ref. |2| due to a reanalysis of 
the data and different choice of q*. Unfortunately, the most direct method, using the heavy- 
light dispersion relation, has large errors. The method using A extracted from heavy-heavy 
mesons is more accurate for aM® < 4.0. For aM® = 7.0 and 10.0, heavy-heavy mesons have 
large discretization errors as these are governed by pa ~ a s Ma, so the corresponding data for 
A are not reliable. To summarize, the best estimate is aM® = 2.31(12) obtained by matching 
M' to the pseudoscalar meson mass, Mb = 5279 MeV. Using M^ in instead of M' gives a 
consistent determination, aM® = 2.21(22), though with larger errors. 

A comparison of the three similar ways of determining M® using the baryon mass 
is presented in Table |5|. Here, and in the following, we use the symbol A^ to represent a 
heavy-light-light A baryon with h labeling the heavy quark. Again, we find that the differ- 
ence between the three methods are significant only for aM® = 7.0 and 10.0. Therefore, we 
determine M® by linearly interpolating the data at the lightest three M® values. The result 
is aM® = 2.5(6) using M^ b = 5624 MeV. This is consistent with the estimate from the meson 
sector; however, since it has much larger errors we do not consider it further. 

The final issue in fixing aM® is related to the fact that our calculation fails to reproduce the 
experimental hyperfine splitting between the B and the B*, as discussed in Sec. [5.4| . Thus, it 
could be argued that determining aM® from the spin-averaged IS* mass (mg+3mB»)/4 = 5313 
MeV should give a better estimate. We find that aM® = 2.32(12) obtained by matching M' to 
the spin-averaged mass is in complete agreement with the value obtained from tub- Henceforth, 
we shall use the value aM® = 2.32(12) for the b quark mass. 
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3 






a E ■ 


aM kin 


aM' 


aMp Crt 


n E ■ 


aM kin 


aM' 




1.6 


0.626(25) 


2.19(31) 


2.28(4) 


2.27(5) 


0.751(14) 


2.54(16) 


2.41(4) 


2.39(4) 


2.0 


0.645(30) 


2.56(40) 


2.66(5) 


2.67(4) 


0.766(16) 


2.93(20) 


2.78(4) 


2.79(3) 


2.7 


0.660(37) 


3.21(53) 


3.35(8) 


3.34(5) 


0.777(18) 


3.62(29) 


3.47(7) 


3.46(3) 


4.0 


0.688(59) 


4.7(11) 


4.68(12) 


4.59(12) 


0.785(28) 


5.02(50) 


4.77(11) 


4.69(9) 


7.0 


0.702(52) 


9.2(11) 




7.44(25) 


0.783(28) 


8.95(56) 




7.52(23) 


10. 


0.726(72) 


16.4(27) 




10.2(4) 


0.783(38) 


14.2(14) 




10.3(3) 



Table 5: E S i m 
Table 0. 



and Ah baryon masses in lattice units. Symbols have the same meaning as in 



4 Mass of the b quark, ^ws^Ws) 

There are two steps needed to determine quark masses from lattice calculations. First, the 
bare quark masses have to be fixed by matching the lattice spectrum to experimental data. 
This has been described in Sec. |3|. Next, one needs to calculate the renormalization constants 
that relate these bare masses to the renormalized mass in the desired continuum scheme. The 
most common scheme is MS and we shall use it here. Standard continuum perturbation 
theory calculations can then be used to convert the result to any other scheme. 

We calculate the MS mass by equating the pole mass on the lattice to that in the contin- 
uum: 

mp i e = Z m Mb = Zconti^mj^ifx) , (2) 
where Z m and Z cont are the lattice and continuum renormalization constants |18j , and \x is the 



scale at which the MS mass is defined. The perturbative series for both Z m and Z cont have 
renormalon ambiguities, therefore so does m po i e . However, in the desired relation, 



m MS^) = ZJLfa) Z m M l 



(3) 



ZcontZm is ambiguity free. 

We calculate m po i e on the lattice in two ways analogous to a previous determination using 
the T system Hl9| . In the first method, we use Eq. (P and write m po i e = M meson — E sim + E 
where M meson is the experimental mass, E sim is measured from the 2-point correlators, and 
E is calculated using perturbation theory. The second method, m po i e = Z m M®, uses the 
perturbative expression for Z m . The quantities Z m , and E , calculated to 0(a s ), are listed in 
Table || for the different values of aM°. The results for the two ways of fixing M 6 ° are given 
in Table ^. 

This pole mass is converted, as in pOU , to Tnj^g(fj) = Z~„ nt {ii)m po i e using continuum per- 



turbation theory for Z~n t (/i) and the Brodsky-Lepage-Mackenzie procedure ]2TJ to set the 
scale for the coupling constant. For /i we choose values between 1/a and n/a, avoiding those 
values where the BLM procedure fails. We then use 2-loop running to get the final result 
m Ms( m Ms)) which, in principle, should not depend on the choice of the intermediate scale 
jj,. These results are also given in Table || where the second error is the spread with respect 



S 



to varying /i, and is indicative of the neglect of the higher order terms in the perturbative 
expressions. 

Our preferred determination of rrijYg(fi) comes from "directly" expanding the product 
Z~J^ t (/j,)Z m in Eq. (|3|) to 0(a s ) [^(J and using the Lepage-Mackenzie procedure |L6] to calculate 



the appropriate scale q* at which to evaluate a s |P0| . The reason for choosing this as the 
preferred method, as explained before, is the cancellation of renormalons in the product and 
the much better value of q*. Continuum (MS) running is then used to convert Tnj^(fi) to 
m lJs^ m 'Ws)- O ur final result, obtained by fixing M° from the spin-averaged M', is 

^ms^ms) = 4.35(10)(;i)(20) GeV, (4) 

where the first error includes statistics and interpolation uncertainty; the second is from the 
uncertainty in the lattice spacing; and the third is the systematic error associated with using 
one- loop perturbation theory We estimate it as being 1 x a 2 s . For typical values of a s , 
depending on the matching scale fi, this is ~ 2.5 — 5%. To be conservative, we assign a 200 
MeV perturbative error to the mass. 

There are two previous lattice determinations of m& using a one-loop matching procedure. 
The NRQCD collaboration [[L^, ^(| |22| has calculated it within the T system, and the APE 



collaboration p3|, |22f evaluates E s i m — Eq for the B meson in the static theory. In addition, the 



APE group has recently extended their matching calculation to two loops ||24j| . These three 
results are 

mf^(mp) = 4.16(5)(15) GeV (NRQCD, 1 - loop), 
= 4.15(5)(20) GeV (APE, 1 - loop), 

= 4.41(5)(10) GeV (APE, 2 - loop). (5) 

While all these results are consistent within errors, a couple of points are in order. First, the 
results of the APE calculation, which is similar to our method 1, suggest that the 2- loop term 
is large. This is consistent with our finding that the aq* for Eq (and Z m ) is small, ~ 0.6. 
Such a small value of aq* is indicative of a large coefficient of the 2-loop term in the bubble 
summation approximation (BSA). Thus in Methods 1 and 2, our estimate of perturbative 
uncertainty in the mass, due to the large value of 1 x a^(q*), is ~ 400 MeV. In our preferred 
direct method, Z~J n± (fj)Z m has no renormalons, and the series is better behaved in the BSA. 
Our estimate of the uncertainty, 200 MeV, is based on the correspondingly larger value of q*. 
To go beyond such an order of magnitude estimate, a two-loop calculation needs to be done 
within NRQCD since the 1-loop calculation shows a strong dependence of the coefficient on 
aM°. Second, we find that the variation of -EWnd = -^sim — -^o with aM° is small, i.e. O(50) 
MeV (see Table ^2|). We estimate that the 0(Aq CD /M) corrections to the APE results are 



of this order. Thus, we expect the systematic error in the APE calculation |24] to be slightly 
smaller than ours. We shall present a more detailed comparison of rrtj^ from the heavy-light 
and heavy-heavy systems on the same configurations in a separate publication B 



5 Heavy- light mesons 

The bare lattice results for meson energies and splittings as a function of k and aM° are 
presented in Tables [7] and |8|. These are first extrapolated/interpolated linearly to Hi and n s , 
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Method to 
fix aM 6 ° 


Pole ma 
Method 1 


ss [GeV] 
Method 2 


ik 

Method 1 


fS mass [G( 
Method 2 


3V] 

Direct 


M' (spin-avg B) 


4.96(1) 


4.97(10) 


4.43(1)(4) 


4.44(10)(4) 


4.35(10)(4) 


M kin (B ) 


4.96(3) 


4.76(41) 


4.43(2)(4) 


4.25(37)(4) 


4.15(38)(4) 



Table 6: Results for the b quark pole and MS masses. Method 1 uses the meson mass and 
E , while method 2 uses Z m and M°. Both methods are described in more detail in [IS]. 



The Direct method is described in the text. The first error quoted is statistical and includes 
interpolation/extrapolation to the physical quark masses; the second is due to the variation 
in the matching scale /i. 

and then to aM 6 ° to obtain estimates for the physical states. (The data are not precise enough 
to include higher order corrections in the fits.) To show the dependence of the mass splittings 
on the heavy quark mass we plot them as a function of 1/M = 4/ (3M#* +Mjj). In this paper, 
we use h to denote a generic heavy quark, H for a heavy-light meson, and an overbar for spin- 
averaged quantities. Where we find a significant M dependence, we quote the intercept (value 
in the static limit) and the slope. In cases where we find no significant slope, we do not show 
the corresponding fits in the figures. In general we find that the slope is ~ Aq CD , i.e. the 
corrections to the static limit are ~ 10% at M&. 

A summary of our results at the b mass is presented in Table |9| and compared with exper- 
imental data in Fig. || We find that the radial and orbital splittings are in agreement with 
the preliminary experimental results. The hyperfine splittings Mb* — Mb and Mb* — Mb s are 
underestimated as will be discussed below. We are able to resolve the P state fine structure 
for the first time on the lattice; previous lattice calculations were done in the static limit and 
found no significant splittings || ^6| . There has been some controversy about the ordering of 



these states in potential model calculations pfl. We find that the Bq is the lightest and B\ is 
the heaviest. Details of the analyses follow. 

In analyzing the mass splittings, we are motivated by the following qualitative picture; the 
mass of a heavy-light hadron is considered to be a sum of: 

• the pole mass of the heavy quark which is ~ 1.5 GeV for the c quark and ~ 5.0 GeV for 
the b; 

• the constituent mass m of the light quarks which is approximately 300 MeV for the u, d 
and 450 MeV for the s quark as inferred from the octet and decuplet light baryons, 
and which we expect to give the biggest contribution to the static binding energy of the 
ground-state hadrons; 

• an excitation energy of the light quark, which, for orbitally and radially excited states, 
we expect to be of the order of Aq C d', 

• the 0(Ag CD /M h ) contributions due to the kinetic energy of the heavy quark and the 
heavy-light hyperfine energy E aH . ai « 45 MeV, inferred from the B* — B splittings; 
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aM° 




a£ sim [lS(%)] 


aE^S^So)] 


aAE(2S - IS) 


a(M H * - M H ) 


aE 


1.6 


0.13690 


0.493(03) 


0.786(23) 


0.293(24) 


0.020(01) 


0.508(03) 


2.0 




0.509(03) 


0.795(23) 


0.286(24) 


0.017(01) 


0.522(03) 


2.7 




0.522(03) 


0.805(24) 


0.283(25) 


0.013(01) 


0.532(03) 


4.0 




0.531(03) 


0.818(24) 


0.287(25) 


0.009(01) 


0.538(03) 


7.0 




0.534(04) 


0.805(25) 


0.271(26) 


0.006(01) 


0.538(04) 


10.0 




0.533(04) 


0.800(26) 


0.267(28) 


0.004(01) 


0.536(04) 


1.6 


0.13750 


0.475(03) 


0.770(27) 


0.295(28) 


0.020(01) 


0.490(03) 


2.0 




0.491(03) 


0.784(26) 


0.293(27) 


0.016(01) 


0.503(03) 


2.7 




0.505(03) 


0.797(29) 


0.292(30) 


0.013(01) 


0.514(03) 


4.0 




0.514(04) 


0.798(28) 


0.284(30) 


0.008(01) 


0.520(04) 


7.0 




0.517(05) 


0.792(28) 


0.275(30) 


0.005(01) 


0.521(05) 


10.0 




0.517(05) 


0.788(30) 


0.272(33) 


0.004(01) 


0.520(05) 


1.6 


0.13808 


0.459(04) 


0.762(33) 


0.303(34) 


0.018(02) 


0.472(04) 


2.0 




0.476(04) 


0.777(35) 


0.301(36) 


0.015(02) 


0.487(05) 


2.7 




0.490(05) 


0.788(31) 


0.298(32) 


0.012(02) 


0.499(05) 


4.0 




0.499(05) 


0.800(38) 


0.302(40) 


0.008(01) 


0.504(05) 


7.0 




0.501(05) 


0.792(34) 


0.292(36) 


0.006(01) 


0.505(05) 


10.0 




0.501(05) 


0.786(35) 


0.285(37) 


0.004(01) 


0.504(05) 



Table 7: aE sim for the IS and 2S mesons is obtained using a two state fit, and a(M H * — Mh) is 
obtained from a fit to the ratio of the correlation functions. The splitting aAE(2S — IS) and 
the spin-averaged energy E = [3E sim (H*) + E sim (H)]/A are calculated within the bootstrap 
process. 



aM° 


K 


uE sim ( 3 P 2 T) 


aAE( 3 P 2 T - 3 P 1 ) 


aAE( 3 P 2 T - 1 P 1 ) 


aAE( 3 P 2 T - 3 P ) 


aAE( 3 P 2 T - 3 P 2 E) 


1.6 
2.0 
2.7 


0.13690 


0.769(08) 
0.774(06) 
0.772(04) 


0.042(11) 
0.042(11) 
0.042(11) 


0.028(07) 
0.028(07) 
0.028(07) 


0.082(11) 
0.078(11) 
0.073(11) 


0.020(14) 
0.020(14) 
0.020(13) 


1.6 
2.0 
2.7 


0.13750 


0.760(09) 
0.765(07) 
0.765(11) 


0.048(13) 
0.048(13) 
0.048(13) 


0.032(09) 
0.032(08) 
0.032(08) 


0.087(12) 
0.083(12) 
0.078(12) 


0.025(16) 
0.025(16) 
0.025(15) 


1.6 
2.0 
2.7 


0.13808 


0.752(10) 
0.757(08) 
0.757(12) 


0.056(16) 
0.055(16) 
0.055(15) 


0.037(10) 
0.036(10) 
0.036(10) 


0.093(14) 
0.088(14) 
0.083(14) 


0.030(20) 
0.029(19) 
0.028(19) 



Table 8: aE s - im and splittings from fits to ratios of correlators for P states. To obtain the 2 + 
P states we spin- average over the 3 P 2 (T) and 3 P 2 (E) states. 
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state (n J F ) 


Lattice 


Expt. 




MeV 


MeV 


heavy-light mesons 


B 1(0-) 
2(0-) 
B* l(l-) 
B* 1(0+) 
B*j 

B* 1(2+) 


5296(04) (jj) 
5895(116) (11°) 
5319(02)(1£) 
5670(37) (1^) 
5770(31)(lf!) 
5822(45) (±H) 


5279 

5860(*) 

5325(1) 

5697(9) 

5779(*)0 
5725-5768(*)|28 


heavy-strange mesons 


B s 1(0-) 


5385(15)(; 6 7 )(^°) 


5369(2) 


2(0-) 


5935(57) (±£)(±g) 




B* s l(l-) 


5412(14)(^)(+J°) 


5416(3) 


i(o + ) 


5742(27) (li 4 )(l^ 5 ) 






5836(25) (11) (i^ 4 ) 


5853(15) 


£; 2 i(2 + ) 


5878(26) (11) (1J 1 ) 





Table 9: Mass estimates in MeV for various meson states. The b quark mass is fixed using the 
spin-averaged B(1S). The first error in the lattice data is statistical (including the statistical 
error in the lattice spacing), the second comes from varying a -1 between 1.8 and 2.0 GeV, 
and for the strange mesons, the third error comes from the uncertainty in the strange quark 
mass. Finite lattice volume effects, which could be large for the excited states, have not been 
addressed in this exploratory study. Preliminary experimental values are denoted by asterisks. 
The lattice results quoted against the Bj and B* sJ states correspond to the spin-average of the 
respective P states, and the experimental numbers are for the unresolved broad resonances. 



Unless stated otherwise, experimental numbers are from the Particle Data Book [22 
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Figure 3: Overview of the B meson spectrum. Circles denote lattice results, dashed lines give 
the range of experimental values [p2]| , and the dotted lines indicate preliminary experimental 
estimates P7|. Errors include statistics and the uncertainty in k s . The variation of a~ l between 



1.8 and 2.0 GeV is not included. 

• and a residual binding energy Ef, e encapsulating the remaining interactions which we 
expect to be small [0(Aq CD /MI)]. 

We accordingly construct different linear combinations of meson and baryon masses to isolate 
individual terms and estimate their size and dependence on the quark masses. 



5.1 B s — Bd splitting 

The spin-averaged splitting between B s and Bd mesons should be dominated by the difference 
of the strange and light quark masses. Our estimate is M Bs — M Bd = 90(9) MeV, 
to be compared to the experimental value 91(3) MeV. The largest uncertainty, the third 
error, comes from setting k s ; estimates using k s (Mk*) are ~ 20% higher, a feature seen in all 
quenched calculations. 

Previous calculations have reported the following results for Mb s — Mb/- 87(1^) (-12) 
MeV H, 86(12)(i£) MeV 0, and 107(13) MeV g| 

(3 = 5.7, 5.9, and 6.1, sees indications of ~ 20% scaling violations between (3 
Averaging the data at the largest two (3 they find 87(7) (4) ( 
87(9)(+|)(±J 9 ) MeV, and the experimental value 90(2) MeV 

In our picture, the heavy quark mass dependence should result from the difference of the 
kinetic and hyperfine energies of the heavy quark in B s and Bd mesons. (In the spin-averaged 
splitting, M b s —M s d i only the difference of the kinetic energies remains.) Therefore, we expect 
this splitting to be independent of the heavy quark mass up to terms of 0((m s — rrid) /M). The 
experimental data show a ~ 10% increase going from the B to the D meson. Our data, given 



The JLQCD calculation |32]] , done at 

5.9 and 6.1. 
1q 9 )- For comparison, our result is 
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in Table |10|, show no significant dependence on the heavy quark mass; however, as shown in 
Fig. |4], they are consistent with the experimental trend. This consistency has also been found 
in Ref. Jl3|, where the heavy quark mass dependence has been studied at higher statistics and 



for a heavy quark mass range between the b and the c. 



AE(H a - H d ) 




lattice units 


MeV 


aM° 


K s (m K ) 


K s (m K *) 


K s {m K ) 


K s (m K *) 


1.6 


0.049(06) 


0.060(09) 


94(10) 


114(14) 


2.0 


0.049(07) 


0.059(10) 


93(13) 


114(17) 


2.7 


0.046(06) 


0.056(08) 


88(11) 


108(14) 


4.0 


0.046(07) 


0.057(09) 


89(12) 


108(15) 


7.0 


0.046(08) 


0.056(11) 


89(15) 


108(18) 


10.0 


0.044(07) 


0.054(09) 


85(14) 


103(16) 



Table 10: Spin- averaged H s 
91(3) MeV. 



Hd splitting as a function of M°. The experimental value is 
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Figure 4: Spin-averaged H s — H^ splitting as a function of the inverse spin-averaged meson 
mass M. The bursts denote the experimental values for B and D mesons. 



5.2 2S - IS splitting 

The raw data for the 2 1 Sq — 1 1 Sq splitting are given in Table 0, and after extrapolation or 
interpolation to Ki and k s , in Table [TTj. This splitting should be dominated by the differ- 
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ence in the kinetic energies of the light and the heavy quarks which give contributions of 
0{.^qcd I m constituent) and 0(Aq CD /M) respectively. With our data, as shown in Table |TT| and 
illustrated in Fig. [|, we cannot resolve any dependence on either the light or the heavy quark 
mass. Our results for B and B s systems are 602(86) (T35) and 559(55) ( jllg) (—12) MeV respec- 
tively, to be compared with the preliminary experimental value, 581 MeV, for the B [27]. In 
the charm sector, the most relevant experimental value is 627 MeV for the D*' ||33|| . 

We do not give results for the spin- averaged splitting 2S — IS, since the signal for the 3 Si 
excited state is less reliable than that for the 1 Sn. 







AE(2 - 1 


L S ) 








lattice units 


MeV 


aM° 


flight 


K s (rn K ) 


K s (m K *) 


flight 


K s {m K ) 


K s (m K *) 


1.6 


0.310(49) 


0.297(27) 


0.294(24) 


593(92) 


570(53) 


564(49) 


2.0 


0.315(54) 


0.294(27) 


0.289(24) 


603(101) 


563(54) 


553(50) 


2.7 


0.313(46) 


0.292(27) 


0.287(25) 


600(91) 


560(55) 


551(53) 


4.0 


0.305(59) 


0.291(30) 


0.287(26) 


585(114) 


558(60) 


551(54) 


7.0 


0.307(53) 


0.280(29) 


0.274(27) 


588(103) 


537(60) 


525(57) 


10.0 


0.299(53) 


0.275(31) 


0.270(29) 


574(103) 


527(64) 


517(62) 



Table 11: 2S — IS splittings extrapolated/interpolated to ki and k s . The preliminary experi- 
mental result for B^ is 581 MeV [f27|, while for B s there is no result as yet. 




0.1 0.2 

l/M [GeV -1 ] 



0.3 



Figure 5: 2S — IS splitting for 1 S states as a function of the inverse spin-averaged meson 
mass. The burst denotes the preliminary 1 5'o experimental value [27 . 
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5.3 IP - IS splitting 



The two main contributions to the spin-averaged IP — IS splitting should be the energy it 
takes to excite the light quark to angular momentum one, 0(Aq C d), and the difference of 
the kinetic energy of the heavy quark in an S-wave and a P-wave light quark background, 
0(Aq CD /M). Our results, shown in Table [12|, are constructed from the raw data given in 
Tables and |. Our estimates are 457(31)(±|f) MeV for the B, and 428(27) (±U)(±H) MeV for 
the B s . 

Experimentally the P states have not been resolved. The P wave resonances Bj (5732) (or 
B**) at 5697(9) MeV and 5^(5850) at 5853(15) MeV are expected to be a superposition of 
the various P states. These are 419 and 484 MeV higher than the corresponding 1 5'o states. 
We use them as estimates of the spin-averaged IP — IS splittings to compare against. 

The variation with either the heavy or the light quark mass is similar to that in the 2S — IS 
splitting. There is a small decrease with increasing light quark mass. The slope, as a function 
of 1/M, is 0.380(202) (ti)(t[j 8 ) GeV 2 for and an almost identical value at shown in 

Fig. § 



AE(P - S) 




lattice units 


MeV 


aM° 


flight 


K s (m K ) 


K s (m K *) 


flight 




K s (m K *) 


1.6 


0.251(13) 


0.238(08) 


0.235(08) 


481(29) 


457(23) 


451(23) 


2.0 


0.244(13) 


0.230(07) 


0.227(07) 


467(27) 


442(21) 


436(23) 


2.7 


0.232(24) 


0.219(08) 


0.216(07) 


446(49) 


420(23) 


414(23) 



Table 12: Spin-averaged P — S splittings. 



5.4 B* - B splitting 

Our results for the hyperfine splitting are shown in Table [13] and plotted in Fig. [7|. A linear 
fit to the B d data gives 0. 138 (38) GeV 2 for the slope and -2(7) MeV for the intercept 
at infinite mass. A zero intercept is consistent with the HQET picture in which the B* — B 
splitting comes from the interaction of the heavy quark spin with the color field, i.e. through 
a a ■ B/(2M) interaction. Our estimates are 24(5) (i 2 ) MeV and 27(3) (I 2 ) (1J) MeV for the B 
and the B s respectively, and AE(B* - B S )/AE(B* - B d ) = 1.19(20)(+|)(±g). These splittings 
are roughly half the experimental values, 46 and 47 MeV respectively. 

An underestimate of hyperfine splittings has also been seen by the previous quenched 



calculations |2j|, [To], |34|, [JT| j32[. The results of the JLQCD calculation [g^] suggest that this is 



not due to scaling violations. Present preliminary unquenched calculations [BJJ do not show 
any significant improvement either, however, the mass of the two flavors of dynamical quarks 
is large, ~ m s . Further work is needed to clarify this issue. 

All hyperfine splittings, including those in the P state and baryon sector, are, to leading 
order, generated by the a ■ B term in the quark action. It has recently been pointed out that 
the coefficient of this term should be larger by a factor of 1.15 — 1.30 |36|, [35|1 . Such a correction 
would bring the quenched results much closer to the experimental values. 
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Figure 6: Spin-averaged P — S splitting of the (circles) and H s mesons (diamonds) as a 
function of the inverse spin-averaged meson mass. The lines denote a linear fit to the data. 



AE(H* - H) 




lattice units 


MeV 


aM {) 


flight 


K s (jn K ) 


K s {m K *) 


flight 


K s {m K ) 


K s {m K *) 


1.6 


0.017(03) 


0.019(01) 


0.020(01) 


32(05) 


37(03) 


38(03) 


2.0 


0.014(02) 


0.016(01) 


0.016(01) 


27(05) 


31(03) 


32(02) 


2.7 


0.011(02) 


0.012(01) 


0.013(01) 


21(05) 


24(03) 


25(02) 


4.0 


0.007(02) 


0.008(01) 


0.008(01) 


13(04) 


16(02) 


16(02) 


7.0 


0.005(02) 


0.006(01) 


0.006(01) 


10(03) 


11(02) 


11(02) 


10.0 


0.004(02) 


0.004(01) 


0.004(01) 


8(04) 


8(02) 


8(02) 



Table 13: H* — H splitting as a function of M°. The experimental results are 45.78(35) MeV 
for B d and 47.0(2.6) MeV for B s . 
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Figure 7: Hyperfine splitting as a function of the inverse spin-averaged meson mass. Circles 
denote the splitting for mesons, diamonds, for H s mesons. For clarity, the diamonds are 
shifted to the right. The burst denotes the experimental value for mesons. The lines are 
a linear fit to the Hd data. 



5.5 P fine structure 

In the jj coupling scheme there are two doublets of P states which are distinguished by the 
angular momentum of the light quark: ji — 1/2 and ji = 3/2. The states in each doublet are 
separated by a spin flip of the heavy quark into a + and a 1 + state for ji = 1/2 {Bq and B*), 
and a 1 + ' and a 2 + state for ji = 3/2 (Bi and B£)- We therefore expect the spin-averages of 
the ji = 3/2 and the ji = 1/2 doublets to be separated by O(Aqcd), and the states within 
each doublet by 0(A 2 QCD /M). 

The experimental situation is as follows. There exists a broad resonance at 5697(9) 
MeV ||22|| , whose spin has not been determined and which is believed to be a superposi- 
tion of various P states. There is also a preliminary experimental result by the DELPHI 
collaboration [27[] for a narrow P state which is 81 MeV heavier than this resonance. Its spin 
is also not resolved, but it is believed to be either J = 1 or J = 2. 

Recently, estimates for individual P states have been obtained by fitting the line shape of 
the broad resonance using phenomenological input based on HQET for the mass splittings, 
decay widths, relative production rates, and branching fractions f28|. Using this method, 
the CDF and ALEPH collaborations obtain a mass of the B\ of ~ 5730 MeV. This result 
seems to be rather insensitive to the assumption about the B\ — B\ splitting, which the 
phenomenological model predicts to be ~ 100 MeV. The L3 collaboration also uses hyperfine 
splittings of 12 MeV as input, but makes no assumption about the splitting between the 
ji = 3/2 and the ji = 1/2 doublet, and obtains slightly higher masses, B\ ~ 5768 and 
B\ ~ 5670 MeV. Our resolution of the P state fine structure is as follows. 
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First, we discuss the + and 2 + states for which the data are shown in Table [14| and Fig. || 
We find that B* -_B* = 155(32) (i? 3 ) MeV and B* 2 - B* = 136(23)(±ig)(±S) MeV - At K h 
the slope versus 1/M is 0.224(70) (^?)GeV 2 and the intercept is 112(33)(ijj) MeV. For B s P 
states, the slope is 0.209(45) (l^)(i^)GeV 2 and the intercept is 97(23) (±|) (±2) MeV. These 
results are a significant improvement over previous values obtained in the static approach, i.e. 
~ 50(100) MeV and ~ 80(75) MeV § for the intercept. 

The situation in model calculations is very unclear. The predictions are model dependent, 
and details like the treatment and the mass of the light quark are significant ||. At this point 
there is no consensus on even the sign of the splitting. 



ae(h; - h*) 




lattice units 


MeV 


aM° 


flight 


K s {m K ) 


K s (m K *) 


flight 


K s {m K ) 


K s (m K *) 


1.6 


0.088(17) 


0.078(11) 


0.075(11) 


168(32) 


149(23) 


144(22) 


2.0 


0.083(17) 


0.073(11) 


0.071(11) 


159(32) 


140(22) 


136(21) 


2.7 


0.078(16) 


0.068(11) 


0.066(10) 


150(32) 


131(22) 


127(21) 



Table 14: if 2 * - H$ splittings. 
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Figure 8: H% — Hq (denoted by octagons) and H* 2 — H* (the diamonds are shifted by 0.01 in 
the x direction for clarity) splittings as a function of the inverse spin-averaged meson mass. 
A linear fit to H 2 — Hq is also shown. 

To study the J = 1 states we used operators with 3 Pi and l P\ quantum numbers in the 
LS coupling scheme as defined in Table |^. The corresponding correlation functions get con- 
tributions from both the physical states. Therefore, at large Euclidean times both correlators 
are dominated by the same lowest state. The masses we extract from short Euclidean time 
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(AEt <C 1) correspond to unmixed states in the LS scheme and are not the physical masses 
||37|| . To get the latter requires a signal in the mixed correlators followed by a diagonalization 
of the 2x2 matrix. Unfortunately, our data does not show a signal in the mixed correlators, 
and therefore we do not have results for the physical J = 1 states. The numbers presented in 
Table |^ under 1 + are those obtained using the 3 Pi correlators. Estimates obtained from the 
1 P\ correlators are almost identical to the center of mass of the 3 P states. 



6 Heavy- light- light baryons 

The heavy- light-light baryons, in the heavy quark limit, can be classified according to the 
angular momentum of the light quarks. At zero orbital angular momentum, the light quarks 
can have total spin si = (anti-symmetric in both spin and flavor) and si = 1 (symmetric 



in both). As summarized in Table 16, there are three states with s/ = 0; udb : usb, and dsb 
which are called the A°, S° and baryons with total spin 1/2. The system with si = 1 
splits up into six hyperfine doublets, each containing states with spin 1/2 and spin 3/2. These 
six doublets are (E+ E£+), (E°,E£°), (Ef^D, (^°,H 6 *°), (E£~,E£~), and (0^") §. The 
pairs of states (E°, A°), and (E' b , E&), do not mix if flavor SU(2) is unbroken. We ensure this in 
our lattice calculation by only analyzing baryons with degenerate combinations of light quarks. 
The raw data are given in Table [T5|. Baryons with a generic heavy quark are denoted as A&, 
2^ etc. To get us and ds combinations we extrapolate linearly in the degenerate light quark 
mass to the average mass (m s + mi)/ 2, which we label K av . A summary of the experimental 
numbers and our lattice results is given in Table [16] and shown in Fig. |9|. 

The UKQCD Collaboration has previously presented a similarly detailed analysis of the 
baryon spectrum 0. They used the tree-level clover action (Csw = 1) at /3 = 6.2 (l/o = 2.9(2) 
GeV) and four heavy k around the charm quark mass. In contrast to our calculation, their b 
spectrum was obtained by extrapolation in 1/M. To facilitate comparison, we summarize their 



results in Table |16|. Within errors these are consistent with our findings, although our results 
are slightly higher and have a slightly smaller light quark mass dependence. An important 
point, as discussed below, is that we are able to resolve hyperfine splittings for the first time. 

The baryon splittings are also analyzed using the phenomenological model discussed in 
Sec. |5|. In heavy-light-light baryons there is an additional light-light hyperfine interaction 
(E av(Tl ), which is expected to be of order Aq C d- 

6.1 A — B splitting 

We first consider the splitting M^ h — (M# + 3M#* ) /4. In this combination, the heavy quark 
mass cancels and there is no contribution from the hyperfine interaction E UH . av Since the 
light quarks are in a ground state with total spin zero, the mass of the extra light quark in the 
baryon gives the dominant contribution. This is borne out by the experimental values: 311(10) 
and 310(2) MeV for the b and c systems respectively, indicating the absence of 0(Aq CD /M) 
contributions from the difference in kinetic energy to the splitting (see Fig. [I]). Our lattice 
data, displayed in Table [T7|, show little dependence on the heavy quark mass, Fig. [LT]. The 
variation with the light quark mass is linear as expected, see Fig. |TT[ Our estimates are 
A b -B = 370(67) (1^) MeV and E b -B s = 392(50) (tj 5 ) MeV. 
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Figure 9: Overview of the b baryon spectrum. Circles denote our lattice results, dashed 
lines give experimental error bounds ||22|| , and dotted lines show preliminary experimental 
results 0,|^. 
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aM° 








0-E sim (E* h ) 


aAE(E* h - E h ) 


1.6 


0.13690 


0.801(10) 


0.852(10) 


0.869(10) 


0.014(02) 


2.0 




0.813(11) 


0.869(11) 


0.883(10) 


0.011(02) 


2.7 




0.823(14) 


0.884(12) 


0.894(10) 


0.008(02) 


4.0 




0.823(20) 


0.899(14) 


0.907(14) 


0.005(01) 


7.0 




0.816(29) 


0.916(20) 


0.921(21) 


0.004(01) 


10.0 




0.806(40) 


0.929(30) 


0.931(32) 


0.003(01) 


1.6 


0.13750 


0.756(12) 


0.822(12) 


0.838(12) 


0.014(02) 


2.0 




0.769(14) 


0.838(13) 


0.851(13) 


0.011(02) 


2.7 




0.781(17) 


0.852(14) 


0.860(11) 


0.008(02) 


4.0 




0.787(27) 


0.866(16) 


0.875(19) 


0.005(02) 


7.0 




0.782(41) 


0.888(23) 


0.895(25) 


0.004(01) 


10.0 




0.776(56) 


0.907(35) 


0.912(38) 


0.003(01) 


1.6 


0.13808 


0.710(17) 


0.795(13) 


0.812(15) 


0.015(03) 


2.0 




0.726(20) 


0.811(14) 


0.825(15) 


0.011(02) 


2.7 




0.739(25) 


0.826(17) 


0.829(13) 


0.008(02) 


4.0 




0.755(41) 


0.835(22) 


0.843(19) 


0.005(02) 


7.0 




0.758(32) 


0.865(29) 


0.876(34) 


0.004(02) 


10.0 




0.766(43) 


0.894(43) 


0.900(45) 


0.002(02) 



Table 15: E sim values for Ah, Eh, and E* h baryons, and E* h — Eh splittings from ratio fits. 



There exist a number of previous results for A& — B, obtained by extrapolating in the 
heavy quark mass, 359(ig)(i^) MeV § and 458(144) (18) MeV ||; in the static limit, 
420(1^°) {tlo) MeV |29|; and with NRQCD on coarse lattices 363(9) MeV (no systematic 
errors quoted). These values are consistent with our result A b — B = 388(68) (^23) MeV. 



6.2 £ - A splitting 

In our picture, the splitting (2E h + 4S^)/6 — depends on E avm , the hyperfme interaction 
between the light quarks, the difference of the binding energies, and of the kinetic energies of 
the heavy quark in each baryon. Experimentally, it is found to be independent of the heavy 
quark mass: (2E c + 4S*)/6-A c = 212 MeV and the preliminary estimate (2S b + 4S^)/6-A fe = 
210 MeV (see also Fig. |l|). These numbers are roughly 2/3 of the Delta-Nucleon splitting (293 
MeV). Such a ratio is obtained in a simple non-relativistic model where these splittings are 
dominated by the light quark hyperfme interaction. The lattice results shown in Table [18] and 
Fig. |l!| are also independent of the heavy quark mass and give 221(71) M eV at M b . 

In the charmed sector the experimental value changes significantly on replacing d with s, 
i.e. (2Hg + 4S*)/6 — S c = 154 MeV. Our lattice results at the b mass also show a decrease 
with 221(71) (ii§) going to 186(51) (±i?)(±? ) MeV, although the difference is not statistically 
significant. 

The UKQCD collaboration || reports E b - A b = 190(±?g) (±|g) MeV and E' b - E b = 
157(lg4) (iii) MeV from extrapolating in the heavy quark mass to the b. Our results for 
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Figure 10: Spin-averaged Ah — H splitting as a function of 1/M. 




Figure 11: Spin- averaged — H splitting as a function of the light quark mass represented 
by the corresponding pseudoscalar meson M%. 
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baryon 


quark content 


experimental 
c ' | b 


1 

b 


Our results 


A-like (s, = 0, j= 1/2) 




(udh) 
(Ish) 


2.285(1) 
2.466 


5.624(9) 


5-64(3) (3) 
5.76(11) (1*) 


5.679(71) 
5.795(53) (±? 5 )(± 5 ) 


E-like (s, = 1, j= 1/2) 


"ft 


(llh) 
(Ish) 
(ssh) 


2.453(1) 
2.574 f2§ 
2.704(4) 


5.797(8) § 


5.77(-)(l|) 
5.90(1^(11) 
5.99(1|) (11) 


5.887(49) ( + ^) 
5.968(39) (11) (If) 
6.048(33) (1*1) (If) 


£*-like (s, = 1, j= 3/2) 


O* 


(llh) 
(Ish) 
(ssh) 


2.519(2) 
2.645 


5.853(8) 


5-78(l & 6 )(lt) 
5.90(1 4 )(1 4 ) 
6.00(11) (1|) 


5.909(47) (1^) 
5.989(39) (til) £T) 
6.069(34) (1^) (If) 



Table 16: Summary of masses in GeV for baryons with quark content shown in column two 
(h denotes a generic heavy quark (c or b), I stands for aword quark). Errors are as explained 
in the caption to Table []. Finite lattice volume effects, which could be large, have not been 
addressed in this exploratory study. Experimental results are given in columns three and 
four. Previous results (UKQCD |J) are in column five. The last column gives results of our 
calculation. 



these splittings are 209(71) and 177(54) (1? ) MeV respectively. 



6.3 £* - £ splitting 

The Y,* h — splitting should depend only on the heavy-light hyperfine interaction E Gh . ai . It is 



therefore expected to be proportional to Our lattice results, shown in Table 19, resolve 

these splittings for the first time. A linear fit to the three lightest M values that bracket M° 
gives — 17(ll)(±i) MeV for the intercept and 0.188(44)(l22)GeV 2 for the slope. However, as 



apparent from Fig. [1^, if the fit is constrained to have zero intercept, then it would have a 



much smaller slope. Based on the assumption that the wavefunction at the origin is similar, 



one expects the slope for the baryon splitting to be 0.75 that for mesons |38[| , which was found 
to be 0.138(38)(lJ7)GeV 2 in Sec. |5T4"| . This expectation does not hold in the charm sector 
where £* - S c « 66 MeV whereas D* - D « 140 MeV. 

The preliminary experimental value is — = 56(8) MeV ||. It is however likely that 
at least one of the states has been misidentified J| , and this number is too large. Scaling the 
experimental value £* — £ c = 66 MeV by M c /M h suggests ~ 20 MeV for this splitting ||, 
We find — £& = 19(7) (!§) MeV; however this could be an underestimate based on the 



general discussion of hyperfine interactions in Sec. pA 

The raw lattice data does not show a dependence on the light quark mass. Experimentally, 
there exists data for strange baryons only in the c sector. The preliminary estimate S* — E' c ~ 
77 MeV is ~ 11 MeV larger than the E* — S c splitting. At the b, heavy quark scaling suggests 
that this difference should be reduced by the factor M c /Mb « 0.3, making it much smaller 
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AE(A h - H) 




lattice units 


MeV 


aM° 


flight 


K s (m K ) 


K s (m K *) 


flight 


K s {m K ) 


K s (m K *) 


1.6 


0.187(25) 


0.263(13) 


0.280(15) 


359(52) 


504(27) 


536(22) 


2.0 


0.191(30) 


0.263(14) 


0.280(15) 


367(60) 


505(31) 


535(25) 


2.7 


0.194(38) 


0.264(18) 


0.280(18) 


372(76) 


506(38) 


536(31) 


4.0 


0.215(59) 


0.266(28) 


0.278(22) 


414(116) 


510(59) 


532(45) 


7.0 


0.228(52) 


0.263(28) 


0.271(30) 


438(104) 


505(60) 


520(60) 


10.0 


0.252(72) 


0.264(38) 


0.267(41) 


484(144) 


507(79) 


513(82) 



Table 17: Splitting between the and the spin-averaged H. The experimental value for 
A h -B d is 310(11) MeV. 



AE(E h -A h ) 





lattice units 


MeV 


aM° 


flight 


K a v(rn K ) 


Kav(jn K *) 


K s {m K ) 


K s {m K *) 


flight 


K av (m K ) 




K s (m K ) 


K s (m K *) 


1.6 
2.0 
2.7 
4.0 
7.0 
10.0 


0.124(29) 
0.119(34) 
0.106(38) 
0.092(59) 
0.126(56) 
0.141(79) 


0.102(21) 
0.099(25) 
0.092(28) 
0.088(44) 
0.118(39) 
0.135(53) 


0.097(19) 
0.095(23) 
0.089(26) 
0.087(40) 
0.116(37) 
0.134(49) 


0.080(14 ) 
0.080(17 ) 
0.078(18 ) 
0.084(29 ) 
0.110(31 ) 
0.129(42 ) 


0.070(14) 
0.071(15) 
0.072(16) 
0.082(23) 
0.106(32) 
0.127(45) 


237(55) 
227(65) 
204(72) 
176(112) 
241(108) 
270(152) 


196(41) 
190(48) 
177(53) 
169(83) 
226(76) 
259(102) 


186(39) 
182(46) 
171(50) 
167(77) 
223(71) 
257(95) 


154(29 ) 
153(33 ) 
150(36 ) 
161(56 ) 
211(59 ) 
248(81 ) 


135(29) 
137(32) 
138(33) 
157(45) 
204(62) 
243(89) 



Table 18: Splitting between the spin-averaged E/, and as a function of M°. K av corresponds 
to setting the light quark mass to (m s +m/)/2. The preliminary experimental value is At = 
210 MeV §. 



AJ5(SJ-E fc ) 




lattice units 


MeV 


aM° 


flight 


K a v(rn K ) 


K av {m K *) 


K s {m K ) 


K s (m K *) 


flight 


Kav(jn K ) 




K s (m K ) 


K s (m K *) 


1.6 


0.016(03) 


0.015(03) 


0.015(03) 


0.014(02 ) 


0.014(02) 


30(07) 


28(05) 


28(05) 


27(04 ) 


26(04) 


2.0 


0.012(03) 


0.012(03) 


0.011(03) 


0.011(02 ) 


0.011(02) 


23(07) 


22(05) 


22(05) 


21(04 ) 


21(04) 


2.7 


0.008(03) 


0.008(02) 


0.008(02) 


0.008(02 ) 


0.008(02) 


16(06) 


15(05) 


15(04) 


15(04 ) 


15(04) 


4.0 


0.005(03) 


0.005(02) 


0.005(02) 


0.005(02 ) 


0.005(02) 


9(05) 


9(04) 


10(04) 


10(03 ) 


10(03) 


7.0 


0.003(02) 


0.004(02) 


0.004(02) 


0.004(02 ) 


0.004(01) 


6(05) 


7(04) 


7(04) 


8(03 ) 


8(03) 


10.0 


0.002(02) 


0.002(02) 


0.002(02) 


0.002(02 ) 


0.003(01) 


3(05) 


4(04) 


4(04) 


5(03 ) 


5(03) 



Table 19: — E^ splitting. The preliminary experimental value for E^ — E& is 56(8) MeV 
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Figure 12: Spin-averaged E^ — A h splitting as a function of the inverse spin-averaged meson 
mass. The bursts denote experimental values for b and c heavy quarks. 



than our resolution. We find E* b - E' b = 19(5) (tf) and Q* b - tt b = 18(4)(±§) MeV. 



7 Heavy- heavy- light baryons 



It is theoretically interesting to study heavy-heavy-light baryons even though it is exceedingly 
hard to produce two overlapping b quarks in experiments. The two heavy quarks are expected 
to bind in a color anti-triplet state whose size is much smaller than Aqqd- It thus interacts 
with the light degrees of freedom to yield a level structure similar to that of heavy-light 
mesons 



41, 42 



In the S'-wave baryons, the total angular momentum of the two heavy quarks is J = or 
1. For identical quarks only J = 1 is possible. There are two different ways to couple the light 
quark spin to this configuration. The J = 3/2 states are denoted as E bb , H^", and f^, and 
the J = 1/2 states as E bb , E bb , and fl bb (the quark content is bbu, bbd, and bbs respectively). 
These are split by a heavy-light hyperfine interaction. Two heavy quarks with different flavor 
can also be in a J = state, and the corresponding baryons are denoted by us as Efy, E' b y, 
and Q'w- The splitting between the spin averaged E&& and the E' w (and the corresponding 
splitting between the Q's) is due to the heavy-heavy spin interaction. This is expected to be 
very small, and to vanish in the infinite mass limit. 



Our raw data are given in Table |20], and the results for E* hh — E hh , extrapolated to m z and 
m s , are listed in Table The data show a strong dependence on the heavy quark mass and 
almost none on the light quark mass. The slope with respect to 1/M is 0.170(42)(±^)GeV 2 
as shown in Fig. |14, and the intercept is — 12(9) (1?) MeV. These results are consistent with 
those for E£ — E^. Both are hyperfine splittings between 5 = 1 diquark and 5* = 1/2 quark 
sub-systems; the difference is whether the S = 1 sub-system is heavy- heavy or light-light. In 
principle the strength of the spin-spin interaction could be different, however, the data suggest 
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Figure 13: T,* h — T, h splitting as a function of 1/M. 

that they are similar. In fact this similarity persists even for h = s where S* — S = 210 MeV 
and £* - £ = 196 MeV. 

If we assume that the spin interaction between the heavy quarks is negligible, then we 
expect {E* hh — ^hh) = l-5(H/j h / — H^). The data shown in Table |2(] indicates a ratio of three 
instead. Our final estimates are 

E bb = 10314(46) (;}?) MeV, tt bb = 10365(40) 0(^ 6 ) MeV, 

~* bb = 10333(55)(; 6 7 ) MeV, Q* bb = 10383(39)Q(^ 2 ) MeV, 

5* 6 - E bb = 20(6)(ijj) MeV, Q* bb - Q bb = 20(4)(^) MeV. 



8 Determination of HQET parameters 

We now present a determination of the HQET parameters A, Ai, and A2. A denotes the binding 
energy of the meson in the limit M° = 00. In the static theory the 0(1/M) corrections to 
this are given by the expectation value of the heavy quark p 2 : 

-\i = ^(B\b(Wfb\B} 7 (6) 

and the expectation value of the chromomagnetic operator: 

\2 = -^(B\ba-Bb\B). (7) 

Thus, to 0(1/M), the relation between the heavy quark pole mass m po i e and the heavy-light 
meson mass is given by: 

M B = rripoie + A + — — — (-Ai + A 2 ) = m voie + E hin d ■ (8) 
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aM° 




nE ■ ( 


Ehh) 


a&E{E* hh - E hh ) 


aAE(E' hh , - E hh ) 


1.6 


0.13690 


0.767( 


08) 


0.015(02) 


0.005(01) 


2.0 




0.788( 


10) 


0.012(02) 


0.004(01) 


2.7 




0.803( 


10) 


0.009(02) 


0.003(01) 


4.0 




0.803( 


14) 


0.007(02) 


0.002(01) 


7.0 




0.767( 


34) 


0.003(02) 


0.001(01) 


10.0 




0.735( 


86) 


0.001(02) 


0.000(01) 


1.6 


0.13750 


0.754( 


10) 


0.015(02) 


0.006(02) 


2.0 




0.777( 


10) 


0.012(02) 


0.004(02) 


2.7 




0.788( 


12) 


0.009(02) 


0.003(02) 


4.0 




0.792( 


15) 


0.007(02) 


0.002(01) 


7.0 




0.746( 


38) 


0.003(02) 


0.001(01) 


10.0 




0.717( 


107) 


0.001(02) 


0.000(01) 


1.6 


0.13808 


0.747( 


12) 


0.015(02) 


0.006(02) 


2.0 




0.768( 


14) 


0.012(02) 


0.004(02) 


2.7 




0.779( 


14) 


0.009(02) 


0.003(02) 


4.0 




0.774( 


19) 


0.006(02) 


0.002(01) 


7.0 




0.733( 


47) 


0.002(02) 


0.001(01) 


10.0 




0.718( 


144) 


-.001(03) 


0.000(01) 



Table 20: E s - im and splittings for heavy-heavy-light baryons. 









^hh ^hh 










lattice units 


MeV 




flight 


K s (m K ) 


K s (m K *) 


flight 


K s (m K ) 


K s (m K *) 


1.6 


0.016(03) 


0.015(02) 


0.015(02) 


31(06) 


29(04) 


29(04) 


2.0 


0.012(03) 


0.012(02) 


0.012(02) 


23(05) 


23(04) 


23(04) 


2.7 


0.009(03) 


0.009(02) 


0.009(02) 


16(05) 


17(04) 


17(03) 


4.0 


0.006(03) 


0.007(02) 


0.007(02) 


11(06) 


13(04) 


13(04) 


7.0 


0.001(03) 


0.002(02) 


0.003(02) 


2(06) 


5(04) 


5(04) 


10.0 


-.002(05) 


0.000(02) 


0.001(02) 


-4(10) 


0(04) 


1(04) 



Table 21: E* hh - E hh splitting. 
In NRQCD one measures -E S im ; from which -Ebind is obtained as 

-EWnd = -^sim — Eq . (9) 

Using the estimates for Eq given in Table ^ -Ebind for the spin- averaged H meson is given in 
Table ^2], and for A/, in Table ^3]. 

We prefer to analyze the dependence of -E^nd on the heavy quark mass in terms of M. The 
reason for this choice is that m po i e is not a physical (measurable) quantity and suffers from 
a renormalon ambiguity. Also, to 0(1/M) the change from m po i e to M is benign, i.e., the 
slope still gives the same Ai and A2 as extracted in conventional HQET analyses. The data 
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Figure 14: E* hh — E hh (circles) and E' hh , — Ehh (diamonds) splittings as a function of the inverse 
spin-averaged meson mass along with linear fits. 



for the binding energy for H, and fits versus M are shown in Fig. [T|. The behavior of A h is 



similar. The results for A and Ai obtained from these fits are also given in Tables U2 and 23 



Note that the slope for spin-averaged cases gives Ai since there is no contribution from the 
chromomagnetic operator. 

Our definition of the parameters A and Ai is perturbative and they inherit a renormalon 
ambiguity from Eq which could be as large as O(Aqcd) in A. With this definition, estimates 
of the HQET parameters are 



A(B) = 375(25) (50) (1^) MeV, 
A(A 6 ) = 895(218) (50) (ijg)MeV, 



-A 1 ( J B)=0.1(3)(l)(^)GeV 2 , 
-A 1 (A 6 ) = -1.7(34)(l)(^)GeV 2 . 



(10) 



We have quoted, as the second error, a systematic uncertainty due to the unknown 0{a 2 s ) error 
in the perturbative expansion of Eq, which we take to be 1 x a 2 s . The third error is due to the 
scale uncertainty. We emphasize that, due to the renormalon ambiguity, these estimates are 
only meant to be indicative and cannot be compared directly with other calculations. 

To remove the uncertainty in A and Ai due to the perturbative estimate of Eq we construct 
differences of binding energies in which Eq drops out. The intercept of a linear fit to the spin- 
averaged Ah — H and — A^ splittings versus 1/M gives 



A(A 6 ) - A(B) 
A(£ 6 )-A(A 6 ) 



415(156)(j^) MeV, 
176(152) (±?„) MeV. 



(11) 



In both cases we find no significant dependence on 1/M. This suggests that the corresponding 
Ai are roughly the same. A similar construction for states with different light quarks gives: 



A(B S 
Xi(B„) 



- A(B d ) 
Xi(B d ) 



*(31)(; 5 d )CT MeV, 
-0.10(28)(^)GeV 2 . 



(12) 
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Ebvad(H) 


aM° 


flight 


K s {m K ) 


K s {m K *) 


1.0 


a no ( 1 o\ 

4Uz(iy j 


AC\CXl '1 /I \ 


oio(oy J 


z.u 


oZO(ZU ) 


410(11 J 


4oy (U 1 ) 


2.7 


378(19) 


466(13) 


485(09) 


4.0 


388(21) 


476(14) 


496(09) 


7.0 


379(21) 


467(13) 


487(10) 


10.0 


375(20) 


460(13) 


478(09) 


A 


375(25) 


458(14) 


477(11) 


-Ai 


0.10(33) 


0.18(14) 


0.20(13) 



Table 22: Binding energies in MeV for the spin-averaged H meson. 



-Ebind(Aft) 




flight 


Kav(jn K ) 


Kav(m K *) 


1.6 


761(60) 


870(47) 


907(40) 


2.0 


692(66) 


798(52) 


833(45) 


2.7 


749(80) 


852(62) 


885(54) 


4.0 


801(120) 


886(94) 


915(82) 


7.0 


817(109) 


888(82) 


912(73) 


10.0 


859(149) 


908(109) 


925(97) 


A 


895(218) 


926(148) 


937(130) 




-1.7(3.4) 


-0.8(2.2) 


-0.5(1.8) 



Table 23: Binding energies in MeV for the A^ baryon. 



Lastly, we estimate A 2 from the slope of the hyperfine splitting calculated in Sec. |5]J. We 
find 

X 2 (B d ) = 0.069(19)(^)GeV 2 , 

X 2 (B S ) = 0.078(12)(^)(± 2 )GeV 2 . (13) 
These parameters have previously been calculated by the Rome collaboration using HQET 



23], ff|. They find 

A(B) = 180(i^)MeV; -\i(B d ) = -0.09(14)GeV 2 ; 

Ai(5 s ) - X 1 (B i ) = -0.09(4)GeV 2 ; X 2 (B d ) = 0.070(15)GeV 2 . (14) 

It is important to note that their definition of A and Ai includes a non-perturbative subtraction 
of the ultra-violet divergence. Thus, the only results that can be compared directly are those 
for \i(B s ) — Xi(Bd) and X 2 (Bd). The experimental values for these two quantities are 

WR , WR x n (M Bs -M B )-(M Ds -M D ) - 
Xi(B.) - X 1 (B i ) = 2 -=——=- = -0.06(2)GeV , 
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Figure 15: -Ebind versus 1/M. 
ME, - Ml 



0.12(l)GeV 2 



(15) 



9 Conclusions 

We have presented an analysis of heavy-light mesons and baryons using a non-relativistic 
formulation (NRQCD) for the bottom quark. Estimates of meson masses with one b quark 
and baryons with one or two b quarks are given in Tables |9| and [16|. Using the B meson 
to fix the b quark mass, we estimate ^jjs^ m Ws) = 

4.35(10) (+i) (20) GeV. This is consistent 



with previous lattice determinations of using the T binding energy [T9|, |2(j, [22| , or HQET 
22, p3l |24]. A more direct comparison will be possible after we extract, using the same set of 



lattices and propagators, m& from the T binding energy. 

A significant feature of our calculation is that we can resolve the P states. We find that 
Mb* < Mb*. Using the interpolating operators based on the LS coupling scheme, we could not 
distinguish between the 1 + and 1 +/ states, as these mix. Also, we resolve b baryon hyperfine 
splittings for the first time on the lattice. 

The mass splittings are analyzed in terms of a qualitative picture based on a non-relativistic 
quark model that is described in Sec. f| We find that the dependence of the splittings on the 
light and heavy quark masses are in agreement with this picture. Quantitatively, the radial 
(25 — 15), orbital (P — 5), S — A, and A — B splittings are found to be within la (~ 20%) 
of the experimental values. 

We are able to resolve hyperfine splittings in both mesons and baryons. The most sig- 
nificant difference from experimental numbers is in the B* — B hyperfine splitting. Such an 
underestimate of hyperfine splittings is a general feature of quenched calculations (light-light, 
heavy-light and heavy-heavy). Another uncertainty associated with the quenched approxima- 
tion is in fixing the strange quark mass. As a result, splittings which are sensitive to the light 
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quark mass have an uncertainty of up to roughly 20% when extrapolated to the strange quark 
mass. 

We have calculated the HQET parameters A, Ai, and A2 for both the B and A and 
Ai have large uncertainties due to the perturbative determination of the shift in the energy 
of the heavy quark, E Q . The differences in these quantities between different hadrons do not 
have this ambiguity and are, therefore, much better determined. 
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